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FUNDAMENTAL GROUP ALGORITHM FOR LOW 
DIMENSIONAL TESSELLATED GW COMPLEXES. 

PIOTR BRENDEL, GRAHAM ELLIS, MATEUSZ JUDA, AND MARIAN MROZEK 


Abstract. We present a detailed description of a fundamental group 
algorithm based on Forman’s combinatorial version of Morse theory. We 
use this algorithm in a classihcation problem of prime knots up to 14 
crossings. 


1. Introduction 

The aim of this paper is twofold. We present the mathematical details of 
the algorithm for computing a presentation of the fundamental group of a 
finite regular CW complex which was introduced in [2] and which is based on 
Forman’s combinatorial version of Morse theory [9]. We also present details 
of the performance of a C++ implementation of the algorithm. To assess 
the performance we systematically apply the algorithm to spaces arising as 
complements of prime knots on fourteen or fewer crossings and use the re¬ 
sulting fundamental group presentation, with standard low-index subgroup 
procedures, to distinguish between the knots. 

1.1. Prior work. The combinatorially defined edge-path group of a con¬ 
nected simplicial complex K, due to Reidemeister, is well-known to be iso¬ 
morphic to the fundamental group 'Ki{K) (see [2Z])- It is also well-known 
that this combinatorial dehnition and isomorphism extends to connected 
regular CW-spaces. In this paper we use the terms edge-path group and 
fundamental group interchangeably as synonyms. Several authors have de¬ 
scribed algorithms for implementing Reidemeister’s edge-path group. Rees 
and Soicher |26j use spanning trees and redundant relator searches in their 
description of an algorithm for hnding a small hnite presentation of the 
edge-path group of a 2-dimensional combinatorial cell complex; they im¬ 
plement their algorithm in GAP [1^ for 2-dimensional simplicial clique 
complexes of graphs. Letscher [21] uses spanning trees and Tietze elimina¬ 
tion/reduction of relators to compute edge-path groups from the 2-skeleta 
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of simplicial complexes arising from knot complements, the knots being pro¬ 
duced from experimental data on protein backbones. Palmieri et al. |25j 
have implemented the edge-path group of simplicial complexes in Sage [2S]; 
the implementation uses the 2-skeleton of the complex and calls gap’s Ti- 
etze reduction/elimination procedures [12j . Kim et al. [TB] describe an al¬ 
gorithm for the fundamental groups of 3-dimensional simplicial complexes; 
their algorithm, which makes use of 3-dimensional cells and the language of 
general CW-spaces, is applied to 3-dimensional tetrahedral meshes arising 
in computer vision. 

We also mention that there is a large literature on the related problem 
of algorithmically determining a collection of shortest generating loops for 
the fundamental group of a space. The case of oriented combinatorial 2- 
manifolds is treated by Ericson and Whittlesey in [10]; their algorithm in¬ 
volves the computation of spanning trees. 


1.2. Our contribution. Recall that a combinatorial vector field on a reg¬ 
ular CW complex X is a partition V of cells of X into singletons and dou¬ 
bletons such that each doubleton consists of a cell and one of its facets, 
that is, one of its faces of codimension 1. The singletons are called critical 
cells, the doubletons are called vectors. The combinatorial vector field is 
acyclic if the facet digraph of X, with direction reversed on edges in V, is 
acyclic. The algorithm we propose (Algorithm |3.1| ) inputs a connected finite 
regular CW complex X, constructs a homotopy equivalent CW complex X' 
with potentially fewer cells, and outputs the finite presentation for vri(X) 
corresponding to the 2-skeleton of X'. The algorithm runs in linear time 
with respect to the number of cells in X (see Theorem 3.2). However, its 
usefulness depends on the size of X'. We have no estimate for the number 
of cells in X' other than it will be no greater than the number of cells in 
X. However, from several numerical experiments we see that X' is usually 
significantly smaller than X. 

At the heart of the algorithm is Forman’s theorem [9l Corollary 3.5] stat¬ 
ing that an acyclic combinatorial vector field V on a regular CW complex X 
leads to another, homotopic, CW complex whose cells are in one-to-one cor¬ 
respondence with the critical cells of V. In order to justify the correctness of 
the algorithm we rework in an algorithmic, recursive spirit the proof of For¬ 
man’s theorem. We actually show that the vectors of V may be quotiented 
out one by one until no vector is left. 

We implemented the algorithm in C-|—|- [Ij. The implementation aims at 
large CW complexes arising from real-world data such as protein backbone 
complements or point-cloud data sets. However, in order to systematically 
study the performance of the algorithm, in this paper we apply it to a family 
of small cubical complexes arising as the complements of prime knots given 
in terms of efficient arc presentations m- We have chosen to work with 
this case, because it provides a large source of examples whose complexity 
can, to some extent, be measured by the number of crossings in the knots. 
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Recall that in [2j we proposed a computable, algebraic invariant of prime 
knots and used it to classify, up to mirror image, all prime knots up to 11 
crossings. The invariant is constructed from low-index subgroups of the fun¬ 
damental group of the complement of the knot with the index not exceeding 
6. The C-|—|- implementation of our algorithm enables fast computation of 
the presentation of the fundamental group. This lets us increase the classifi¬ 
cation in [2] from 801 knots up to 11 crossings to 59937 prime knots up to 14 
crossings in the tabulation provided by Hoste, Thistlethwaite and Weeks in 
m and available from m- The C-|—|- implementation shifted the computa¬ 
tional barrier from the presentation of the fundamental group computations 
to the algebraic invariant computations. The latter are performed in GAP 
[I2] and HAP [7]. 

The organization of the paper is as follows. In Section we gather nota¬ 
tion and preliminaries. In Section we present a recursive approach to the 
discrete Morse theory. In Section]^ we study various reduction techniques. 
In Section we present the algorithms. In Section we apply the results of 
the paper to knot classihcation. 


2. Preliminaries 


2.1. Notation. Let N := {1, 2, 3, ...} stand for the set of natural numbers, 
M for the set of reals and let I denote the interval [0,1]. We write / ; X-e^Y 
for a partial map from X to Y, that is a map defined on a subset dom f C X, 
called the domain of /, and such that the set of values of /, denoted im/, 
is contained in Y. Given a topological space X we write cl A for the closure 
of A C A. 

For n G N we write 


5" 


Bn 

:= (x G M” \\x 

n—1 

:= (x G M” X 

n—1 

:= {x=(xi,X2,. 

n—1 

:= {x=(xi,X2,. 

o 

Bn 

:= (x G M” X 


^ 72—1 

1 


a:|| < 1}. 


Xn ^ 0 }) 

Xn ^ 0 } 5 


Note that in particular = [—1,1] and 5° = { — 1,1}. We extend this 
notation to n = 0 by setting 


■= rO := {*}, 

S-^ := 0, 


where * is a unique point not belonging to any other set considered in this 
paper. 

Let A be any set. An oriented element of A is := (a, e) for e G { — 1,1} 
and a G A. We denote the set of oriented elements of A by A. We identify 
a^ with a and we say that a and a~^ have mutually inverse orientation. 
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By a path in a topological space X we mean a continuous function 0 : 

= [—1,1] —)• X. We identify the oriented path 9^ with 6 and 9~^ with 
the path 

e{-t) G X 

We refer to 9~^ as the inverse path of 0. 

2.2. Directed multigraphs. Recall that a direeted multigraph is a quadru¬ 
ple G = {V, E,e-i,e+i), where e_i,e+i : E ^ V are given maps. The set 
V is the set of vertiees of G, the set E is the set of edges of G and the maps 
e_i, e+i assign to each edge its initial and terminal vertex respectively. 

A path in G is a sequence of oriented edges p := such that 

Geki^k) = (efc_|_i) for A; = 1,... n — 1. The initial vertex of the path p 

is e_ej(ei) and the terminal vertex of the path rj is ee„(e„). The path r/ is a 
cycle if the terminal and initial vertices coincide. The inverse of the path rj 
is 77-1 := 

Let C be a path in G with the same initial and terminal vertices as an 
edge z G E. Moreover, assume that neither z nor z~^ appear in Let r] be 
a path in G. A (z, C)substitution in rj, denoted Sz^(^{rj), is the path obtained 
from rj by replacing every occurrence of 2: in 77 by C and every occurrence of 
z~^ in 77 by 

2.3. Quotient spaces. We recall from [6j some facts concerning quotient 
topological spaces. Given an equivalence relation R on a topological space 
X, we denote by the equivalence class of x G X and hy X/R the set of 
all equivalence classes of R. Let kr : X 3 x is [x]h G X/R be the canonical 
projection. We drop the subscript R whenever the relation R is clear from 
context. The quotient topology on X/R consists of subsets A C X/R such 
that k~^{A) is open in X. 

Proposition 2 . 1 . [6l Proposition 2.4.14] Assume R is an equivalence re¬ 
lation on X and S is an equivalence relation on X/R. Then, X/R/S is 
homeomorphic to X/Rs, where Rs is an equivalence relation on X defined 
by xRsy if and only if ks{kr{x)) = Ks{KR{y)). 

Given a partial continuous map / : Y sX such that dom / is closed 
in Y we define X Uf Y as the quotient space {X \JY)/Rf of the disjoint 
union X U T by the equivalence relation Rf which identifies x G dom / with 
f{x) G X. For z G X U Y we denote by [z]f the equivalence class of the 
relation Rf. Note that if / is empty, then X DfY coincides with X UY. 
In the special case when Y = B^ and dom / = we say that X UjY 

results from X by gluing an n dimensional cell via the attaching map f. 

The relation R is called closed if the map kr is closed. The following 
proposition is a special case of m Proposition 2.4.9]. 

Proposition 2 . 2 . The relation R is closed if and only if for any open U in 
X the set 


Ur:={xGU\ [x]rCU} 
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is open in X. 


□ 


The following proposition is an easy consequence of Proposition 2.2 


Proposition 2.3. Assume Y <Z X is closed and R is a closed equivalence 
relation in Y. Then R := i2Uidx is a closed equivalence relation on X. □ 


Proposition 2.4. Assume X is a Hausdorff topological space and R is a 
closed equivalence relation in X with compact equivalence classes. Then 
X/R is a Hausdorff space. 

Proof: Let [x] / [y] be two equivalence classes of R. Since they are 
compact and X is a Hausdorff space, we can find disjoint open sets U, V in 
X such that [x] C U and [y] C V. By Proposition [2^ and Vr are open in 
X. It follows that k{Ur), k{Vr) are disjoint, open in X/R. Obviously, they 
are neighbourhoods respectively of [x] and [y]. Hence, X/R is a Hausdorff 
space. □ 

Given H C X we denote hy X/A the quotient space X/Ra, where 

Ra '■= { (x, y) G X"^ I X = y or X G H and y G H. } 

2.4. Finite CW complexes. In this paper by a CW complex we always 
mean a finite CW complex. Here we recall its definition. Let iL be a Haus¬ 
dorff topological space. A finite CW structure on K is a pair {IC,{(pa}a£K.) 
such that /C is a finite family of subsets of K and: 

O 

Each element cr G /C is a subset of K homeomorphic to H"" for some 
n > 0. The subset a is referred to as an n-cell. The number n is 
called the dimension of a and denoted |it|. The dimensions provide 
a filtration := U { ^ ^ ^ I ^ J }■ 

The family /C is a decomposition of K, i.e. K = [JJC and any two 
different elements of JC are disjoint. 

O 

Each ip(j is a continuous map ip^ : —)■ K mapping homeo- 

morphically onto a and C . The map p>(j is called 

the characteristic map of a. 

A finite CW complex (cf. O Section 1.3]) is a triple (AT,/C, {(^o-jo-eA:) such 
that (/C, {pa/a&ic) is a finite CW structure on K. It is convenient to slightly 
abuse the language and refer to K as the CW complex, assuming that the 
CW structure is given implicitly. 

The CW complex is called regular if each pa ■ —)• AT is a homeomor- 

phism onto its image. 

By the boundary of a cell a we mean the set Bd a := cla\a. We say that 
a cell T is a face of a cell u if r n cl cr 7 ^: 0 . The face is proper ii t a. The 
face is a facet if |r| = |cj| — 1 . A 0-dimensional cell is called a vertex. 

We refer to the collection of facets of a cell cr as its combinatorial bound¬ 
ary and we denote it by bdcr. It is also convenient to consider the collec¬ 
tion of cells whose facet is cr. This collection is called the combinatorial 
coboundary of a and is denoted cbdcr. Note that if ct is a 1-cell, then 


(CWO) 

(CWl) 

(CW2) 
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bdcT = (/?(!)}. In the case of a 1-cell these two facets (or one if 

<^(— 1 ) = (/?(!)) are referred to as the endpoints of a. 

A cell is top-dimensional if it is not a proper face of a higher dimen¬ 
sional cell. The CW structure and CW complex are called pure if all top¬ 
dimensional cells have the same dimension. We will refer to the closure of a 
top-dimensional cell as a toplex. It is easily seen that a CW complex is the 
union of its toplexes. 

A union L := (J £ of a subfamily C <Z JC such that L is closed in K is called 
a subcomplex of K. A subcomplex is a CW complex with its CW structure 
consisting of C with the same characteristic maps as in JC. Note that, in 
particular, every vertex is a subcomplex. Also, (J is a subcomplex, called 
the j-sJieleton of K. 

Since by m Theorem 1.4.10] the closure of a cell in a regular CW complex 
is a subcomplex, we have the following proposition. 


Proposition 2.5. If cells cr,T in a regular CW complex satisfy r n cl cr 7 ^ 0 , 
then r C cl(T, i.e. t is a face of a. □ 


As a consequence of Proposition 2.5 we obtain the following proposition. 


Proposition 2.6. Assume t is a proper face of a in a regular CW complex 
X. Then t is a facet of a face of a. 


Proof: Set n := |(t|. Since X is regular, Bd cr is homeomorphic to 
Let It := |cr| — |r|. Since r is a proper face of a we have A: > 0. We 
proceed by induction on /c. If A; = 1 the claim is obvious. Thus, assume 
A: > 1. Let x € t and consider a shrinking sequence of open balls in Bdcr 

O 

centered at x. All these balls are homeomorphic to and none of these 

balls is contained in r, because otherwise |r| = n — 1 and A: = 1. Thus, 
there exists a sequence {xj} C Bd cr \ r converging to x. By passing to a 
subsequence, if necessary, we may assume that the sequence is contained in 
a cell a' C Bdcr. It follows that r n clu' 7 ^ 0 and, by the regularity of X, r 
is a face of a'. The conclusion now follows from the induction assumption, 
because Ict'I — |r| <Jt. □ 


Proposition 2.7. A |in[ Proposition 1.5.17] j A finite CW complex is a metriz- 
able topological space. □ 


Note that if h : —>■ is a homeomorphism, then h maps homeo- 

O O 

morphically B"' onto B^ and onto S^~^. In consequence, if is a 

characteristic map, then so is (pa o h. Thus, we can freely substitute home- 
omorphisms into characteristic maps without loosing the CW structure. In 
particular, in the characteristic map pa we can replace by any set home¬ 
omorphic to B'^. 

Let fj be a top-dimensional cell of K. Then K\a \s easily seen to be a 
subcomplex of K and K results from K\a hy gluing via the attaching 
map := 
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Let a = be an admissible ordering of the cells of /C, i.e. an ordering 

of all cells in /C such that each cell is preceded by its faces. For j = 1,2,... k 
set 

j 

T,j := (Tj. 
i=l 

It is not difficult to prove that each Sj is closed in K, hence a subcomplex 
of K and Sj results from by gluing the cell aj via the attaching map 

One easily verifies that an admissible ordering always exists on a finite 
CW complex. Moreover, when L is a subcomplex of K, then one can choose 
the ordering in such a way that L = S; for some I < k. Thus, we have the 
following standard result. 

Proposition 2.8. Each finite CW-complex can be constructed recursively 
by attaching a cell to a subcomplex. □ 

We will use also the following fact which is a special case of |20l Theorem 

11 .11]. 

Theorem 2.9. Assume Xi, X 2 are homotopy equivalent spaces and h : 
Xi X 2 is a homotopy equivalence. If maps fi : —)• Xi are such 

that f 2 = h o fi, then the spaces Xi and X 2 U /2 B^ are homotopy 

equivalent. □ 

2.5. Oriented CW complexes. Our fundamental group algorithm of a 
CW complex X is based on a theorem reducing the computation to a group 
generated by the 1-cells with relators provided by the 2-cells. To present 
the theorem in detail, we recall the definition of the orientation of an n-cell 
a in X. For n = 0 it is just an element of {—1,1}. For n > 0 it is an 
equivalence class in the homotopy equivalence relation of the characteristic 
map ip a- : —)• (clcr, Bdu) considered as a map of pairs. By [T3l 

Corollary 2.5.2] there are precisely two orientations for each cell. An oriented 
CW complex is a CW complex with a fixed orientation for every its cell. 

Assume K is an oriented 2-dimensional CW complex. By an edge in K we 
mean an oriented 1-cell. An edge r, as each 1-cell, has two endpoints (^r(—1) 
and Pt{1)- The orientation allows us to distinguish between them. The first 
is called the initial vertex of r and denoted e_(r), the other is called the 
terminal vertex of r and denoted e+(r). Thus, the CW structure of the 1- 
skeleton of K may be considered as a directed multigraph. A characteristic 
map of an edge is a characteristic map pr of the associated 1 -cell which 
induces the chosen orientation of the edge. 

Consider now a 2-cell a. There is a closed path ri, r 2 ,..., in the directed 
multigraph of edges such that the homotopy class of the attaching map 
6 * 0 - contains a map 0:5^—)- such that = Uj=i h 
characteristic map of the edge Tj (see m Proposition 3.1.1 and 3.1.2]). The 
loop ti,T 2 , ... ,Tn is Called the homotopical boundary of a. We have the 
following theorem. 
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Theorem 2.10. / |13[ Proposition 3.1.7 and Theorem 3.1.8]j Assume K is 
a connected CW complex with precisely one vertex. Then, the fundamental 
group of K depends only on the 2-skeleton of K. Moreover, up to an iso¬ 
morphism it is the group generated by the edges of K with arbitrarily selected 
orientation of K and homotopical boundaries of all 2-cells as relators. 


3. Discrete Morse theory 


3.1. Reduction pairs. We say that a facet oq of an n-cell ai is regular if 
ip ax, the characteristic map of oi, satisfies 

(i) ipax{S'^~^) is a subcomplex of K, 

(ii) iPax\B”- is the characteristic map of ao, 

(hi) (faxiBI) nao = 0. 

Note that in a regular CW complex each facet is automatically regular. 

We say that a := {oq, ai) is a reduction pair in a CW complex K if oq is 
a regular facet of oi. The following proposition is straightforward. 

Proposition 3.1. If L is a subcomplex of K and a is a reduction pair in 
L, then a. is a reduction pair in K. □ 

3.2. Discrete vector fields. A discrete vector field V on the CW structure 
/C is a partition of K. into singletons and doubletons such that each doubleton, 
when ordered by dimension, forms a reduction pair. By a critical cell of V 
we mean a cell a € )C such that {a} is a singleton in V. We denote the set 
of critical cells of V by crit V. By pairing the lower dimensional cell with the 
higher dimensional cell in every doubleton of V we may and will consider V 
as an injective partial self-map V : tC-e^tC such that JC is the disjoint union 
of domV, imV and critV and for each a G domV the pair {a,V{a)) is a 
reduction pair. We will call such a reduction pair a vector of V. 

With V we associate a directed graph, Gy. Its vertices are the cells of 
/C and edges go from each cell a to each of its facets r with orientation 
reversed if u = V(r). If Gy is acyclic, the discrete vector held V is called 
acyclic. In this case the graph Gy induces a partial order on the cells in JC. 
In particular, for each cell a £ 1C there is a minimal cell less than or equal 
to (T in this partial order. We say about such a cell that it is inferior to a. 

Assume now that 77 is a regular CW complex with a CW structure 1C. 
An algorithm constructing an acyclic discrete vector held on 1C is presented 
in Table It is a simplihed version of the algorithm proposed in [MS], 
based on the method of coreductions [23] . 

For a CW structure 1C let M(/C) denote the maximal cardinality of bdci 
and cbdcj over all u G /C. 


Theorem 3.2. Assume Algorithm 3.1 in Tahle^is called with 1C containing 
a CW structure of a non-empty, regular, connected finite CW complex K. 
The algorithm always stops, returning an acyclic vector field on 1C. More- 
on return C contains the respective set of critical cells with precisely 


over, 


one critical cell in dimension zero. If the data structure used to store the 
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Algorithm 3.1. discreteVectorField(CW structure /C) 

V — L — C := 0 ; 
while /C 7 ^ 0 do 
if L = 0 then 

n := min { g | /C^ / 0 }; 
move an element r from /C” to C; 
for each n G (/C n cbdr) \ L do enqueue u in L; 
else 

cj:=dequeue(L); 

if /C n bdfj = {r} then 

remove r and a from 1C and insert (r, a) ioV] 
for each u G (/C n cbdr) \ L do enqueue u in L; 
else if /C n bd O' = 0 then 

for each u G (/C n cbd cr)\ L do enqueue u in L; 

return V; 


Table 1. Discrete Vector Field by coreduction method 


CW complex provides iteration over the boundary and coboundary of a cell as 
well as removing a cell in 0{M) time with M := M(/C), then the algorithm 
runs in 0 (M^ card/C) time. 

Proof: Observe that if an element t £ 1C contributes its coboundary 
elements to L as an element of a coreduction pair (r, ct), then it is removed 
from 1C, so it may never again contribute in this role. Similarly, if an element 
a with zero boundary contributes its coboundary elements to L, it may never 
again do so, because it is not in the coboundary of any element, so it may not 
reappear in L. Therefore, each element of S may appear in L at most 2M 
times. It follows that the while loop is passed at most 2Mcard/C times. In 
particular, the algorithm always stops and its complexity is 0 (M^ card/C). 

Let ICj, Cj and Vj denote respectively the contents of variable 1C, C and 
V on leaving the jth pass of the while loop for j > 0 and the initial values 
of these variables for j = 0. We will show by induction in j that Vj is an 
acyclic vector field on K, with ICj U Cj as the set of its critical cells. This is 
obvious for j = 0, because Vb = 0, so every cell is critical. Thus, assume the 
claim holds for k < j. It is obvious that Vj is a discrete vector field on 1C. 
We only need to prove that Vj is acyclic. This is evident if Vj = Vj-i, thus 
assume Vj / fj-i. Then Vj = Vj-i U {(r, cr)} for some reduction pair (r, a). 
Assume Vj is not acyclic. Then, there exists a cycle so,to, si,ti, ... Sn,tn, sq 
in Gvj such that G Vj and ti is a facet of S(j+i) modn- particular, 

tn G bdsQ. Since by induction assumption is acyclic, precisely one of 
the pairs {ti, Si) is the pair appearing in Vj on the jth pass of the while 
loop. Without loss of generality we may assume that this is the pair {to, sq). 
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Thus, we have JCj-i n bd so = {^o}- In particular, to G Since the pair 

(ti,si) G it must have been included in the vector field on an earlier 

pass of the while loop, say kth. pass with k < j. Thus, IC^-i Hbdsi = {ti}. 
But to G bdsi, therefore to 0 K,k-i- However, the sequence (/C*) is clearly 
decreasing, hence to 0 a contradiction. This proves that Vj is acyclic 

for every j. In particular, the vector held returned by the algorithm is 
acyclic. Since after the last pass of the while loop /C = 0, its critical cells 
are all stored in C. 

We still need to prove that C contains precisely one vertex. First observe 
that for each j the union of the cells in set Aj := Cj U domV^- U imV^- 
is a closed subcomplex of K. Indeed, this is obviously true for j = 0, 
because = 0- Now, arguing by induction, we see that Aj+i is either 
empty, or a critical cell in Cj+i whose boundary is in Aj, or a doubleton 
{r. O'} such that (r, a) is a reduction pair satisfying K.j n bdci = {r}, i.e. 
bdcj\^j = {r}. Thus, in each case \JAj is closed, because [_}Aj-i is closed 
by induction argument. 

Observe that on the very hrst pass of the while loop a vertex, say v, 
is removed form /C and moved to C. Thus, we only need to show that this 
never happens again. Assume the contrary. Let w be another vertex moved 
to C and assume this happens on the jth. pass of the while loop. Since 
K is connected, there is a path v = vq, eo,vi, ei, ... em,Vm = vj such that 
each Si is a 1-cell with endpoints in vertices Vi and Uj+i. Without loss of 
generality we may assume that Vi G C implies Vi = v or Vi = w. Since 
\_]Aj-i is closed and Vm 0 we see that also We will 

show that also Vm-i 0 H iiot) then Cm would have been placed on 

L and then removed together with Vm as a reduction pair before the jth 
pass of the while loop. Arguing by induction, we see that vo 0 Aj-i, a 
contradiction. □ 

3.3. Quotients of reduction pairs. For a reduction pair a we write Kq. := 
cloi, (fa ■■= (fai, |a| := |ai|, Ka ■= K+ := k\ := 

K\ai \ ao- Set 

K-a ■= { T G/C I r C iLo }, 

/Cq, . A-q \ oi}. 

Lemma 3.3. We have Ka = = IJ^a- particular, Ka 

and K^ are subcomplexes of K. 

Proof: Obviously \JJCa C Ka- To prove the opposite inclusion take 
an X G Ka- If x G oi, then x G IJ^Q) because oi C JCa- Otherwise, 
X G By assumption (i) of the definition of a regular facet we 

have X G r for some cell r C But C /Cq, hence 

X G U JCa - This proves the hrst equality. 

We have 

O 

Ka = U U Mb'+) = ai U ao U iL+. 
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But ai n = 0 by the definition of a cell and ao n = 0 by the 
definition of the reduction pair. Therefore, 


K'^ = Ka\ai\ao = 1J/C+, 

which shows that Ka is a subcomplex of /C. □ 

Letx = (yi,X2):BH^i?H-ix / be a homeomorphism such that 

x{b'^-) = X X I, 

x(i?t') = i?H-ix{0}. 

Set Xa := 

Define the map /„ : Ka —)• K'^ by 


fa{x) 


X 


if X G «! U ao, 
otherwise. 


It is not difficult to check that fa is well defined and continuous. 


Proposition 3.4. The map fa is a retraction of Ka onto K^, i.e. fa\K+ ~ 
id^+. Moreover, it is a homotopy inverse of the inclusion map K'^ C Ka- 


Proof: The fact that fa is a retraction is an elementary calculation. To 
see that fa is a homotopy equivalence, define a homotopy h : Ka x I —)> Kf^ 
by 


h{x, t) 


Xa{XliTaHx)),tX2{Tj{x))) 

X 


if X G ai U ao, 
otherwise. 


It is not difficult to check that h is well defined, continuous and ho = /a, 
hi = id^+. □ 

OL 

Define an equivalence relation Ra in K by 


Ra := { (x, y) G I X = y or X, y G Ka and fa{x) = fa{y) }■ 


Let [x\a denote the equivalence class of x with respect to Ra and let Ka '■ 
iL 9 X I— >■ [x]a G K/Ra denote the quotient map. 


Proposition 3.5. The space K/Ra is Hausdorff. 


Proof: Note that [x]q = /“^(/q(x)), therefore the equivalence classes 
of R are compact. Since K is Hausdorff and Ka is closed in K, in view of 
Proposition |2.4| and Proposition |2 . 3| it is enough to verify that the restriction 
of Ra to Ka is closed. By Proposition 2.2 it suffices to show that for any U 
open in Ka the set Ur is open in Ka- Assume this is not true. Then, for some 
X G Ur no neighbourhood P of x satisfies V C Ur- Thus, Proposition |2.7| 
allows us to select a sequence x^ —> x such that [xn]_R U- Let yn G [xn] \D. 
Passing to a subsequence, if necessary, we may assume that yn ^ y & K\U- 
Since yn G [xn], we have faixn) = faiVn)- Thus, /q,(x) = faiv)- It follows 
that [x]ij = [y]R- In consequence y £ U, a contradiction. □ 


Set Kit := n 


'a\K 


+ . 

oc 
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Lemma 3.6. The restriction k\ := : Ka — ?■ K/Ra is a continuous 

bijection. 

Proof: The map k\ is continuous as a restriction of a continuous map. 
To show that k\ is surjective, take an x £ K. If x 0 Ka, then Kq,(x) = 
[x\a- Hence, assume x £ Ka- Since fa is a retraction onto we have 
fa{x) = faifaix)). ThuS [x]a = [/a(x)]a = K+(/q,(x)), which shoWS that 
K'f is surjective. Assume in turn that [xi]^ = [x 2 ]a for some xi,X 2 G K. 
Then /q,(xi) = /q,(x 2 ). If either xi or X 2 is not an element of Ka, then 
both classes are singletons, hence xi = X 2 - Thus, assume that xi,X 2 G K'^. 
Since fa is a retraction onto K^, we have xi = fa{xi) = faix 2 ) = X 2 , which 
shows that nf is injective. □ 


Lemma 3.7. The restrictions and Ka\K\Ka homeomorphisms 

onto their images. 


Proof: By Lemma 3.6 both maps are continuous bijections. Since Ka is 
compact and K/Ra is Hausdorff, the first restriction is a homeomorphism 
onto Kf^. To see that the other restriction is a homeomorphism observe that 
equivalence classes in iL \ Ka are singletons, therefore i<^^{Ka{U)) = U for 
any U C K\Ka. In particular, for any open set U C K\Ka its image Ka{U) 
is open, which shows that also the other restriction is a homeomorphism. □ 


3.4. Collapses. We call the quotient space K/Ra an a-collapse of K and 
write briefly Ka for the quotient map kr^. 

By performing a collapse we do not change the homotopy type and stay 
in the class of CW complexes, as the following theorem shows. 


Theorem 3.8. The quotient space K/Ra is a CW-complex with CW struc¬ 
ture ({NQ(cr)}^6>c\{«i,ao}) ({«a o V?-T}aeA{“i,ao})- Moreover, the complexes 
K and K/Ra are homotopy equivalent. 


Proof: First recall that by Proposition |3.5| K/Rg is a Hausdorff space. 
Let a £ 1C \ {ai,ao}- In order to prove property (CWO) of the dehnition 
of CW complex it suffices to show that Kaicr) is homeomorphic to a. By 
Lemma 3.3 the set K/^ is a subcomplex, therefore either cr 
K\K, 


C iPa or (T C 

Thus, by Lemma |3.7[ in both cases KQ-(cr) is homeomorphic to a, 
as required. Properties (CWl) and (CW2) are straightforward. 

By Proposition |2.8| complex K may be obtained from Kg by consecutively 
gluing a sequence {aj) of cells via the attaching maps 6^.. It is not diffi¬ 
cult to observe that applying to Kg/Rg the same sequence of gluings but 
with attaching maps Kg o Og. we obtain the complex K/Rg. Therefore, the 
conclusion follows by induction from Theorem |2.9| and Theorem |3.8[ □ 


Proposition 3.9. If L is a subcomplex of K, then Ka{L) is a subcomplex 
of K/Rg. □ 
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Let V be an acyclic vector field on JC. Choose ao,a:i G V such that 
ai = V(ao)- Then a := (ao,«i) is a reduction pair. Define a partial map 
V« : K/R^ ^ K/R^ by 

domVa := { Kq((t) | cr G dom V \ {ao} }, 

Va(Ka(fT)) := Ka{V{a)). 


Theorem 3.10. If V is acyclic, then Va is a well defined discrete vector 
field on K/Ra and it is also acyclic. 


Proof: To show that Va is a well defined discrete vector field on K/Ra 
consider /3o G domV \ {oq} and let fii := V{fio). Then, /3 := (/?i,/3o) is 


a reduc tion pair in K 
Lemma 


3.6 


Moreover, /3 / a. Thus, /3i,/3o £ Since by 
Ka is injective on Ka, we see that Va is injective and domVo n 
imVo = 0 . Hence, we only need to show that Kaifio) is a regular facet 
of Kgifi i) in JCa- Note that (i) follows from Proposition 3.9 and (ii) from 
Lemma 


3.6 


_ To prove (hi) assume the contrary. Then, n 

0- Hence, we can choose an x G and a y £ fio such that x y 


and Ka{x) = Ka{y). Since Kg is bijective on Kg, we have x,y £ Kg and 
either x or y is not in Kfi. But y £ uqU ai implies fio = uq or fio = ai and 
contradicts a fi. Thus, x 0 i.e. either x G ai or x G ao- Consider 
X G «! . Then ai C 0 which implies Kg C K'^ and 


0 fion K+ cfion K+ cfion Kg c fion iL+, 


a contradiction. Consider x G ao- Then, /3o n and ao H K'^, i.e. ao is a 
facet of fii and fio is a facet of ai. Thus, ao, ai,/3o, ,0i, ao is a cycle in Gy, 
a contradiction again. This completes the proof that Vg is a well defined 
discrete vector field on K/Rg. 

We still need to prove that Vg is acyclic. Assume it is not. Then, we have 
a cycle 

Ka(/3o), K«(/3i), Ka(/3f), • • • 


such that /3* := (/3Q,/3j) 7 ^ a, Ka(/3Q)) is a regular facet of Ha{j3\) and 
Kq,(/3q'''^)) is a facet of Ka{l3\). Since I3q,I3\ C Kpi and K^i is a compact 

subset of k\, we see from Lemma 3.6 that Kg\x i is a homeomorphism, 


hence /3 q is a regular facet of j3\. If cl/3j n Kg = 0 for all i, then 

( 1 ) 


is a cycle in Gy, because Kg restricted to iP \ Kg is a homeomorphism onto 
its image. Hence, assume that is not a facet of for some i but 
C cl Ka(fil). This is only possible if ao is a facet of fi\ and is a 
facet of ai- Thus, we can insert ao, ai between and /3j in 0 obtaining 
a cycle in Gy, a contradiction. □ 

A sequence a = (a^, a^,..., a™') of reduction pairs in JC is called ad¬ 
missible if the associated collection {a^, a^,..., a™'} of reduction pairs is a 
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discrete vector field. Obviously, given a discrete vector field V, any ordering 
of its elements is admissible. Assuming V is acyclic and fixing an order¬ 
ing a, we can use Theorem 3.10 to recursively apply collapses of reduction 
pairs in a. By Theorem 3.8, the resulting space is a CW complex homotopy 
equivalent to the original complex K. We call it the Morse complex of the 
vector field V and denote K/V. This is some abuse of language, because 
the construction depends on the ordering a but we assume the ordering is 
implicit. Summarizing, we get the following fundamental result of discrete 
Morse theory. 


Corollary 3.11. Let K he a finite CW complex and V an acyclic discrete 
vector field on the CW structure K, of K. Then, the complexes K and K/V 
are homotopy equivalent. 


□ 


4. Reductions. 

4.1. External collapses. We say that the reduction pair a is free if is 
a sub complex of K. 

Theorem 4.1. If a is free then the complexes K and k\ are homotopy 
equivalent. 


Proof: We know from Lemma 


3.6 


that the restriction k is a con- 

tinuous bijection onto K/Ra. Thus, it is a homeomorphism, because 
as a subcomplex, is compact. Hence, the conclusion follows from Theo- 

□ 


rem 


When a is free, we say that k\ is an external collapse of K. Otherwise 
we say that K/Ra is an internal collapse of K. On the algorithmic side 
external collapses are advantageous, because to compute them it is enough 
to remove the reduction pair from the list /C of cells of K. Thus, it is 
convenient to make as many external collapses as possible. We say that K 
collapses onto the subcomplex L if there exists a sequence of subcomplexes 
L = Lq C Ti C ... C Lra = K such that Lj-i is an external collapse of Lj. 
Note that for convenience we admit the case when m = 0, i.e. L = K. 

An immediate consequence of Proposition is the following proposition. 


Proposition 4.2. If K collapses onto a subcomplex L then the inclusion 
L <Z K is a homotopy equivalence. □ 


4.2. Collapsible subcomplexes. Complex K is called collapsible if it col¬ 
lapses onto a vertex. 

Theorem 4.3. Assume L is a collapsible suhcomplex of K. Then, K and 
K/L are homotopy equivalent. 

Proof: Let Lm C Lm-i C ... C Lq = L be a sequence of subcomplexes 
of L such that Lj is an external collapse of Lj_i by a reduction pair 
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and Lm = {*} consists only of a vertex. Let kj : Lj —>■ Lj /R^j be the 
quotient map and let k* := k™ o o ■ ■ ■ k} . By Proposition |3 .1 1 each 
is a reduction pair in K. Thus, by Theorem |4.1| K is homotopy equivalent 
to 


iL/i?„i/i?Q,2 ■ ■ ■ /Ra 


and by Proposition 2.1 this space is homeomorphic to K/R^, where R^ is 
a relation on X defined by xR^y if and only if k*(x) = K^{y). But k* is 
constant on L and a homeomorphism on K\L. Hence, K/R^ and KjL are 
homeomorphic, and consequently K and K/L are homotopy equivalent. □ 


Theorem 4.4. Assume K is a regular CW complex and L is a subcomplex 
with CW structure C C 1C. Then, K collapses onto L if and only if there 
exists an acyclic discrete vector field on 1C such that crit V = C. 


Proof: Assume first K collapses onto L. Then, we have subcomplexes 
L = Lq <Z Li <Z ... C. Lm = K such that Lj-i is an external collapse of 
Lj resulting from removing a reduction pair . It is straightforward to 
verify that the collectionja^, ... a”} is an acyclic discrete vector field on 
/C whose set of critical cells is precisely L. 

To prove the opposite implication assume V is an acyclic discrete vector 
field on K, such that crit V = C. Let n be the cardinality of /C\£. We proceed 
by induction on n. For n = 0 the conclusion is trivial. Thus, consider the 
case n > 0. Then, we can choose a non-critical cell in cr G /C. Since V is 
acyclic, we can choose a cell r which is inferior to a, that is r is a minimal 
cell less than or equal to a in the partial order induced by V. We claim that 
r is not a critical cell of V. To see this consider a path r = tq, ri, ..., = cr 

in Gy- Since r is inferior to a, such a path obviously exists. If r is critical, 
then we consider the smallest index k such that Tk-i is critical and Tk is not 
critical. Since the partial order induced by Gy goes against the facet partial 
order only for the vectors of V, we have Tk G domV and ta, C clr C L, 
because r, as a critical cell is a subset of L and L as a subcomplex is closed. 
It follows that Tk is critical, a contradiction. Therefore r G domV. Let 
K' := K \ t\ V(r). To prove that K' is closed, assume the contrary. Then, 
there exists an x G cliL' \ iL'. Hence, x G r n V(r) and x G cl cr' for some 


cell o' C K'. By Proposition 2.5 r is a face of o' and by Proposition 2.6 
r is a facet of a face o" of o'. Thus o" precedes r in the partial order 
induced by Gy, and consequently r is not inferior to o. This contradiction 
shows that K' is a subcomplex of K and consequently K collapses onto K'. 
By induction assumption K' collapses onto L. Therefore, K collapses onto 
L. □ 

As an immediate consequence of Theorem 4.4 we obtain the following 
corollary. 


Corollary 4.5. A regular CW complex is collapsible if and only if it admits 
a discrete vector field whose set of critical cells consists of precisely one 
vertex. □ 
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4.3. Redundant configurations. Let iL be a regular CW complex, L its 
subcomplex and K, and £ the corresponding CW structures. We refer to 
the pair (iL, L) as a configuration. We say that the configuration (iL, L) is 
redundant if there exists a discrete vector field V on /C such that crit V = £. 
We say that subcomplexes L, M of iL form a redundant decomposition of K 
\{ K = L\J M and the configuration (M, L n M) is redundant. 

Theorem 4.6. Assume subcomplexes L and M of a regular CW complex K 
form a redundant decomposition. Then, the inclusion L C K is a homotopy 
equivalence. 


Proof: Let V be a discrete vector field on M such that crit V consists of 
the cells in LC M. Extend it to all cells in K by making all cells not in the 
domain of V critical. It follows from Theorem 4.4 that K collapses onto L, 
thus we get the assertion from Proposition 4.2 □ 


The advantage of Theorem 4.6 is twofold. We may use it to reduce the 
complex without changing its homotopy type by removing a subcomplex 
from a redundant decomposition. We may also simplify the complex by 
quotienting out a maximal collapsible subcomplex. We can construct such 
a subcomplex L C K hy starting from a vertex and recursively adding to 
L a subcomplex M such that L and M form a redundant decomposition of 
L U M. However, to make it work in practice we need a quick method to 
test for redundant configurations. 

This is possible in the setting of tessellations discussed next. 


4.4. Tesselated spaces and lattice complexes. We say that a CW struc¬ 
ture /C on a pure, regular CW complex X, finite or not, is a tessellation if 
the closures of any two top-dimensional cells cti, ct 2 in JC are isomorphic, i.e. 
there exists a homeomorphism h : cl ai —?■ c1iT 2 mapping cells onto cells. A 
CW complex is a tesselated CW complex if its CW structure is a tessellation. 
Obviously, every pure subcomplex of a tesselated space is a tesselated space. 

A special way to define a tessellation is to consider a lattice L C M"', i.e. an 
additive subgroup of generated by some choice of n linearly independent 
vectors Vl = {ui, • • • , Vn}. Any v € L determines a Dirichlet-Voronoi cell 

Dl{v) := {x G I ||u — x|| < ||rc — x|| Vrc G £}. 

For X G MT set 

xUx) := {v £ L \ X G Dl{v) } 

and define an equivalence relation R C M” x M” by xRy if and only if 
Xl{x) = Xiiv)- One can show that the equivalence classes of this relation 
provide a CW structure of M”" which is a tessellation. We call it a lattice 
tessellation and denote it JCl. When Vl consists of n orthogonal vectors, the 
sets Dl{v) are re-dimensional hypercubes and the tessellation is the classical 
cubical decomposition of M”". By identifying M” with the hyperplane 

{(xi,... ,x„+i) G : xi H-h x^+i = 0} 
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and taking as Vl the set vi = (1 — re, 1,1 ,..., 1,1), n 2 = (1,1 — R, 1, • • •, 1,1), 

..., Vn = (1,1,1,..., 1 — re, 1) we obtain a tessellation of M” consisting of 
re-dimensional permutahedra. 

A pure, finite CW complex whose CW structure is a subfamily of a lattice 
tessellation will be called a lattice CW complex. Lattice CW complexes may 
be easily and efficiently stored in bitmaps. 

Fix a lattice L. Two top-dimensional cells va. tCi are called neighbours if 
the intersection of their closures is non-empty. Consider a top-dimensional 
cell a G ICl and a lattice CW complex K with its CW structure JC <Z JCl- 

Observe that K\a and cl a are also lattice CW complexes. Their inter¬ 
section is a subcomplex of X. We will call it the contact complex of a and 
K. Note that this makes sense both if cj C iL and a Ci K = 0. 

Let Nk{(t) denote the set of neighbours of cr contained in K. Then 

(iL\cj)ncliT = |^{cl(Tnclr| T G Nk io'),T / a}, 

i.e. the contact space is uniquely determined by the neighbours of a in K. 

The toplex cl a is said to be redundant in K if the configuration (cl a, K\ 
a n c\a) is redundant. If a has k neighbours, then there are 2^ configura¬ 
tions to test for redundancy. Since all toplexes in a lattice tessellation are 
isomorphic, tests done for one toplex may be reused for other toplexes. Even 
better, for reasonably small number of neighbours we can test all configura¬ 
tions of a toplex for redundancy and store them in a lookup table for later 
quick redundancy tests. The number of neighbours for a cubical tessellation 
in M" is 3” — 1. The number of neighbours for a permutahedral tessellation 
in M”' is — 2. Thus, the lookup tables may be computed and stored in 
less than 1GB of memory for cubical tessellations up to dimension 3, requir¬ 
ing 2^® bits or 2^^ bytes in dimension 3 and for permutahedral tessellations 
up to dimension 4, requiring 2^*^ bits or 2^^ bytes in dimension 4. 


5. Reduction algorithms. 


5.1. Shaving and collapsing. Two reduction algorithms based on redun¬ 
dancy tests are presented in Tables and The first, often called shaving, 
removes redundant toplexes as long as they exist. The second constructs a 
possibly large collapsible subcomplex. 

Theorem 5.1. Algorithm \5.1\ applied to a collection of toplexes representing 
a lattice space always stops and returns a homotopy equivalent subcomplex 
with no redundant toplexes. 


Proof: The conclusion follows immediately from Theorem 4.6 


□ 


Theorem 5.2. Algorithm \5.S\ applied to a collection of toplexes represent¬ 
ing a lattice space always stops and returns a collapsible subcomplex. The 
resulting quotient complex is homotopy equivalent to the complex on input. 


Proof: The conclusion follows immediately from Theorem 4.6 and The¬ 
orem 14.31 □ 
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Algorithm 5.1. shaving(collectioii of toplexes K) 

Q:=0; 

while there exists a redundant toplex a in K do 
remove a from K] 
enqueue Nk{(t) in Q; 
while Q ^ 0 do 
(T=dequeue(( 5 ); 
if fj is a redundant in K do 
remove a from K; 
enqueue Nxify) in Q; 

return K; 


Table 2. Shaving 


Algorithm 5.2. collapsibleSubset(collection of toplexes K) 
C ■=Q := 0 ; 

choose any toplex cr in A and insert it into C; 
enqueue ™ Q; 

while Q 7 ^ 0 do 
(T:=dequeue(( 5 ); 
if a is redundant in C do 
insert a into C; 
enqueue Nx\(j{a) in Q; 
return C; 


Table 3. Collapsible subset construction 


5.2. C-structures. In order to obtain an algorithm based on acyclic vector 
helds and Theorem 2.10| it is convenient to consider the following concept. 

Let C be a finite set decomposing into the disjoint union C = Cq U 
Cl U 6 * 2 . By a C-structure we mean a quadruple (C, e_i, e+i, d) such that 
(Co, Cl, e_i, e+i) is a directed multigraph and d : C 2 —)• Cycles(Ci), where 
Cycles(Ci) stands for the set of cycles in Ci. 

We write |c| = j if and only if c G Cj. 

For c G C we define its boundary by 


0 if |c| = 0 , 

{e_i(c),e+i(c)} if |c| = 1 , 

{ ei, 62 ,... e„ } if |c| = 2 and d(c) = ■ ■ ■ e^. 


bdc:= < 
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With every C-structure C we associate a 2-dimensional CW complex K 

O 

defined as follows. For each c G C we have a formal cell (c) ^{c}- 

However, to simplify notation, in the sequel we identify (c) with c. Let 

:=[j{c^C\ |c| < j}. 

For c G Cl and t G we define the characteristic map (pc'- ^ by 


:= J iftGHi, 
[et(c) otherwise. 


To define the characteristic map pc '■ B^ 
d{c) = take an x G Set 


for a c G C 2 assume 






if X GH^, 

if X = and —k < t < —(k + 1). 

n — — n N / 


It is not difficult to verify that the cells c together with the characteristic 
maps Pc define a CW complex. We denote it by CW(C). Conversely, given a 
finite, regular, 2-dimensional CW complex K it is straightforward to obtain 
a C-structure C such that CW(C) is isomorphic to K. 

Given a reduction pair (ao,Q;i) in CW(C) we define a C-structure C as 


follows. Let C' := C \ {ao, cii}. To define e' : C( —)■ Cg and d! : C2 


Cycles(Cj), consider first the case |ai| = 1. By the definition of a reduction 
pair do is a regular facet of ai. This implies that bd ai consists of precisely 
two elements one of which is ao- Let be the other. Set 


e'eic) 

d'{c) 


I Og if ec(c) = ao, 

[ee(c) otherwise, 
d(c). 


Consider the case |ai| = 2. Since ao is a regular facet of oi, the facet oq 
appears precisely once in the cycle d(ai) = ■ ■ ■ e'pp. Without loss of 

generality we may assume that ao = ei. Let ag := Then 

ag and ao have the same initial and terminal vertex. Moreover, ao does not 
appear in ag. Set 

d\c) := S’„o,„ 2 (d(c)), 

where denotes the ao,a5 substitution in d{c). 

It is not difficult to verify that (C', e'_g, e'_)_g, d') is a C-structure. We call 
it the a-collapse of C. One easily verifies the following proposition. 

Proposition 5.3. The CW complex CW(C') associated with the a-collapse 
C is homeomorphic to the quotient complex CW(C)/i2a. □ 


Algorithm |5.3| computing a presentation of the homotopy group of a reg¬ 
ular CW complex K represented as a collection of top-dimensional cells is 
presented in Table 
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Algorithm 5.3. fundGroup(collection of toplexes K) 
K := shaving(A); 

A := collapsibleSubset(A); 

C := C-structure of K/A] 

V := discreteVectorField(C'); 
for each a G 1/ do 

assign to C the a-collapse of C; 
return (Ci,d(C 2 )); 


Table 4. Fundamental group presentation of a regular CW complex. 


Theorem 5.4. Assume Algorithm 15.31 is called with a collection of top cells 
of a CW structure of a lattice space K. It always stops and returns a pre¬ 
sentation, up to an isomorphism, of the fundamental group of K. 


Proof: The theorem follows immediately from Theorem 5.1, Theorem 5.2 


Theorem|3.2[ Theorem|3.8[ Proposition|5.3[ Theorem 3.10 and Theorem 2.10 


□ 


6. Application to knot classification. 

6.1. Knot classification. Recall that two knots K, K' C are equivalent 
[Sam] if there exists a homeomorphism / : —)• such that K' = f{K). 
They are orientation equivalent |22j or isotopic ffH if the homeomorphism 
is an ambient isotopy of knots, that is it is isotopic to the identity. Gordon 
and Luecke m Theorem 1] prove that K, K' are equivalent if and only if 
the complements \ iF and \ K' are homeomorphic. Thus, topological 
invariants of the complements of knots may be used to distinguish between 
non-equivalent knots. They cannot be used to distinguish between non¬ 
isotopic knots. But, by a result of Fisher [ 8 ] two equivalent but non-isotopic 
knots differ only by the mirror image. 

Let T be an invariant of topological spaces, that is T{X) = T(Y) if X and 
Y are homeomorphic. Then, obviously T{M.^\K) = r(M^\iF') if K and K' 
are equivalent. We say that T classifies a family of knots K, if for any two 
non-equivalent knots K,K' € JC we have T(M^\iF) 7 ^ T{M.^\K'). Note that if 
the family JC is finite and the invariant T is algorithmically computable, then, 
at least theoretically, one can verify if T classifies 1C in finite computation. 
A natural candidate for a good T is the fundamental group. Unfortunately, 
by the results of Novikov [21] and Boone [T], in general it is not possible 
to decide algorithmically from finite presentations of two groups if they are 
non-isomorphic. However, this finite presentation may be used to compute 
some further, algorithmically decidable, invariants. This is t he path we used 
in [2]. Recall that given a knot iF C M^, using Algorithm|5.3|we can compute 
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c 

3 

4 

5 

6 

7 

8 

9 

N{c) 

1 

1 

2 

3 

7 

21 

49 

c 

10 

11 

12 

13 

14 

15 

16 

N{c) 

165 

552 

2176 

9988 

46972 

253293 

1388705 


Table 5. The growth of the number of prime knots N{c) 
with the number of crossings c. 


a presentation of the fundamental group Gk ■= tt{D^ \ Mk), where is a 
cubical ball in containing the knot K and Mk is a cubical neighbourhood 
of K containing K as a deformation retract. Using this presentation and 
the GAP software library, we can compute some algebraic invariants of Gk- 
In particular, given integers n,c > 1, m > 0 we study in [2] the invariant 
It is defined by 

j[n,cM^K) := {H^iShc+iS,!.) : S<Gk,\Gk : < n}, 

where Hm{G,'L) denotes the integral homology of group G and is de¬ 
fined recursively by 71 G := G, 7 fc+iG := [ 7 fcG, G]. Trying to balance 
the computational expenses between the cost of computing the invariant 
(AT) and its ability to distinguish between knots we chose the in¬ 
variant /"'(AT) := l[”’^’^l(Ar) with c = m = 1 and varying n. Using a GAP 
implementation of Algorithm |5.3| we provide in [ 2 ] a computer assisted proof 
of the following result. 

Theorem 6.1. The invariant I^{K) distinguishes, up to mirror image, be¬ 
tween ambient isotopy classes of all prime knots that admit planar diagrams 
with eleven or fewer crossings. □ 


6.2. Classifying index. Using a G-|—|- implementation of Algorithm |5.3| 
and an improved strategy of computing the classifying invariant we are able 
to strengthen Theorem 6.1 (see Theorem 6.3 below). The computational 
challenge of this extension lies, in particular, in the exponential growth of the 
number N{c) of prime knots with c crossings. Prime knots up to 16 crossing 
were tabulated by Hoste, Thistlethwaite and Weeks in [16]. The number of 
such knots for a given n is presented in Table In our computations we 
used the Hoste-Thistlethwaite tabulation available from m- 

There are only 801 prime knots up to 11 crossings but 59937 prime knots 
up to 14 crossings. The additional challenge is in the cost of computing the 
invariant 7”, particularly when n grows. An observation we made is that for 
most knots relatively low values of n suffice to distinguish them from other 
knots. But, some knots require relatively high values of n. Thus, in order 
to decrease the computational cost, it makes sense to take as the classifying 
invariant the pair (n,/"(AT)) for a possibly low n selected individually for 
each knot. Assuming the invariant classifies a family of knots 1C, given a 
knot K £ 1C we refer to n in the invariant (n, I'^{K)) as the classifying index 
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Algorithm 6.1. classify(collection of knots fC) 

Q := empty queue; //to process (A, n) 

A := empty set; / / {K, n) added to Q 

I := empty dictionary; // I[i] is a list of {K,n), s.t. i = I^{K) 
C := empty dictionary; // C[K] final invariant for K 
D := empty dictionary; // D[i\ = A iff C[K\ = i 
A := 0; // classifying index 
foreach A G /C do 
enqueue {K, 2) in Q] 
insert (K, 2) into A; 
while Q is non-empty do 
{K, n):=dequeue((5); 

N := max(A, n); 
i ■.= r{K)- 

if i is not a key in I then 
I[i\ := empt list; 
append {K,n) to list I[i]-, 

D[i] := A; 

C[K] := (n,i); 
else 

append {K,n) to list /[i]; 
remove key i from D; 
foreach {K, n) G I[i\ do 
(*) remove K from C; 

if {K, n + 1) ^ A then 
(**) enqueue (A, n -|- 1) in Q; 

insert {K, n -|- 1) into A; 
return (C,D,N); 


Table 6. Classification algorithm. 


of K. By the classifying index of the family K, we mean the maximum of 
all classifying indexes of individual knots K € 1C. The idea of classifying 
indexes leads to Algorithm |6.1| presented in Table Note that the while 
loop in Algorithm 6.1 is potentially infinite. But, we have the following 
theorem. 


Theorem 6.2. Assume Algorithm \ 6.1\ starts with a collection of knots K, 
on input and stops. Then, the triple {C, D, N) on output, consisting of 
dictionaries C, D and an integer N, has the following features. 

(a) The integer N is the classifying index of JC. 

(b) The keys of C are all knots in fC. 
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(c) The keys of D constitute a eollection of invariants I^{K), one per 
eaeh knot K £ fC. Given an invariant i in the keys of D, the as¬ 
sociated value D[i] is the unique knot K in JC such that C[K] = 
{n,I'^{K)) = {n,i). Then, the number n is the classifying index of 

K. 

(d) Both dictionaries C and D are injective. In particular, for any two 
knots K,L £JC the invariants C[K] and C[L] are different. 

Proof: Given a variable x in Algorithm 6.1 we let x^^'l denote the state 


of X after I iterations of the while loop. In particular, for I = 0 this is the 
state just before entering the while loop. We will show by induction that 
for every / > 0 we have the following properties: 

(ii) the sequence { is not descending, 

(hi) / {K'd ')for each V < I, 

(iv) for each K £ JC the value of C^^'l[K] is set if and only if AT G JC\C^^\ 
where 

£« ■.= {K£JC\ Gg«}, 

(v) C(') and are injective, 

(vi) the value D^^^i] is the unique K £ JC such that 

(2) C^^\K] = {n,r{K)) = i. 

First observe that dictionaries are empty, is equal to 

because of the first for loop, contains only pairs {K, 2) for every 
K £ JC, = JC, is not defined, = 2 and > 2 for j > 1. In 
particular, properties (i) (vi)| are trivially satisfied for I = 0. Thus, fix f > 0 
and assume the conditions are satisfied for k < 1. 

To see |(i)| observe that in course of the /th iteration of the while loop 
we add an element to A if and only if we add it to Q and we never remove 
elements from A. This proves [(I)| 

By th e ind uction assumption for k < 1. Hence, in order 

to prove (ii) we only need to show that The inequality is 

obvious when = 2. Thus, assume that > 2. Then, 

and consequently also enter the queue Q in pair with a respective knot 
inside the while loop. Hence, there exist k and k' such that I > k > k', 
+ 1 and ^ + 1. From the induction assumption 

^ therefore = jii^) -|- 1 > + 1 = 

To prove (hi) let V < I he such that This 

obviously cannot happen if ^ = 2. Thus, for I there exists a, k < I and for 
f there exists a k' < f such that 


and 


(i^(n,„(n)Gg(^')\gT). 
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Without loss of generality we may assume that k' < k. In the fcth iteration 
line (**) is executed with {K^^\n ^ However, 


a contradiction. _ 

Now we prove (iv) Assume to the contrary that there exists a K G 
such that C^^'^[K] is not set. Then, on the Zth pass of the while loop line {*) 
is executed with = K. Indeed, ii K ^ then is set and 

the algorithm executes line (*) by the induction assumption. If AT G 
then K = and the algorithm executes the else branch, because 
is not set. In particular, (AT, n) with n = is processed as one of the 
elements of in the second foreach loop. If the algorithm executes 

line (**), then K G )C^^\ a contradiction. Otherwise (AT, n + 1) G A^^ and 
consequently (A', n + 1) G for some I' < 1. Because of the assumption 
(AT, n + 1) also is not in . Hence, K = AT^^ ^ and n + 1 = ^ = n 


for some I" such that I' < I" <1. This contradicts (ii) 

Assume in turn that C^^\K] is set and K G If is set, then 

by the induction assumption K 0 Thus, the algorithm executes line 

(**) in the /th iteration with ATi = K. In particular, it also executes line 
(*), a contradiction. If is not set, then AT^^) = K. Let n := vM\ 

Since K G and the algorithm stops, there is a A: > / such that = 
K. Without loss of generality we may assu me th at k is minimal with this 
property. By (ii) we have and by (iii) > rS-^\ It cannot be 

jii^) = j7,(0 _)_ because on the Ah pass of the while loop the else branch 
is not executed. Thus, + 1. It follows that there exists a p < k 


such that = K and + 1. Therefore, and by (ii) 

p > 1. Consequently, the number k, contrary to its choice, is not minimal. 

It follows easily from the induction assumption that 


Now we prove (v) 


if is not injective, then there exists a k < I such that 


It means that 

(n('),i(')) = (n(^),iW). 

We set ^[Ar] only in the then branch of the first if statement. Hence, 
^(0 ^ But in the kth iteration of the while loop we set C'[A'(^)]. 

Therefore, G We do not delete elements from I. Thus, since 

k < I — 1, we have G a contradiction. 

Again from the induction assumptions, if is not injective, then there 
exists a k < I such that 

a:(o = D(0[j(0] = = K^k)^ 

and Since we have z^^^ = = /"(AT^^^) = 

a contradiction. 
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c 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

N{c) 

2 

2 

3 

3 

3 

3 

5 

5 

6 

6 

7 

7 


Table 7. Classifying indexes N{c) for the family HTW<c 
with c < 14. 


To prove (vi) let i, K be such that = K. Since D is injective, we see 

that K is unique. If then Q follows from the induction 

assumption. Otherwise, the algorithm sets D in the Ith iteration. Thus, it 
also sets C[K] to (n,i). 

To finish the proof first observe that assertion (a) is obvious. Let I* de note 
the final pass of the while loop. Then, = 0. Thus, we get from (iv) 


that the keys of C are all knots in JC. This proves assertion (b). Assertion 
(c) follows from (vi) and assertion (d) from |(v)] □ 

Given a family JC of prime knots and an integer c we denote by JCc the sub¬ 
family of these knots in JC whose minimal planar diagram requires c crossings 
and we set JC<c '■= Ui<c^*- HTW denote the family of 1 701 936 prime 
knots tabulated by Hoste, Thistlethwaite and Weeks in |16] and available 


from |19j . We have the following extension of Theorem 6.1 


Theorem 6.3. The invariant I^[K) distinguishes, up to mirror image, be¬ 
tween ambient isotopy classes of all prime knots in HTW<i 4 and the clas¬ 
sifying index of this family is 7. 


Proof: Apply Algorithm |6.1| to the family HTW<i 4 . The algorithm stops, 
returning 7 as the classifying index of the family. Hence the result follows 
from Theorem [Q □ 

The same way we can determine the classifying index for HTW<c with 
c < 14. The results are presented in Tablej^ Figure [^presents the grouping 
of classifying indexes of individual knots. Tables and indicate that 
for most knots their classifying index is relatively low. 


6.3. Limitations of the method. Recall that during the computations 
of the invariant of a knot K we first build a finite presentation of the 
fundamental group G of \ AT and construct the family S'^{G) of the 
subgroups of index n oi G. 

Since each subgroup of index n of a finitely generated group G corresponds 
to a homomorphism G ^ Sn, the number of such subgroups is bounded from 
above by (n!)® where g is the number of generators in G. We cannot bound 
the number from below. Thus, we can influence the computational time 
of constructing S'"'(G) only via g. Therefore, since g is not greater than 
the number of critical cells of the constructed discrete vector fields, it is 
important to minimize this number. The problem of computing S^{G) is 
well studied in the literature. A good source of references is [3]. 
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Figure 1 . The image shows knots separation in HTW<i4 by 
the invariant /”■. The middle white circle represents all knots 
in HTW<io. The first ring counting from inside presents arcs 
grouping knots with the same invariant with classifying index 
2 . The length of the coloured arc represents the group size. 
The following rings show groups for the classifying indexes 3 , 
4 and 5 . We can see how the groups split. Also, it is visible 
that the classifying index exceeds 3 only for a small number 
of knots. 


The problem of computing the optimal number of critical cells is NP- 
complete HZ! In the construction of the discrete vector field we use greedy 
algorithms. They usually behave very well for homology computation, where 
the difference of one more or less critical cell has no significant consequence 
for the total computational complexity. However, in the case of the P in¬ 
variant, one more critical cell may increase the total computation time from 
minutes to days. The most time-consuming knot we encountered is K14,33437 
presented in Figure]^ For this knot the computations of P{K) vary between 
a few minutes and 52 days, depending on the method used to construct the 
discrete vector field (coreductions, reductions) and the order of cells in the 
data structure. The best method and order varies from knot to knot. Since 
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c\n 

2 

3 

4 

5 

6 

3 

1 ( 100 %) 

1 ( 100 %) 

0 ( 0 %) 

0 ( 0 %) 

0 ( 0 %) 

4 

1 ( 100 %) 

1 ( 100 %) 

1 ( 100 %) 

0 ( 0 %) 

0 ( 0 %) 

5 

2 ( 100 %) 

2 ( 100 %) 

2 ( 100 %) 

2 ( 100 %) 

0 ( 0 %) 

6 

3 ( 100 %) 

3 ( 100 %) 

3 ( 100 %) 

2 ( 67 %) 

0 ( 0 %) 

7 

7 ( 100 %) 

7 ( 100 %) 

6 ( 86 %) 

3 ( 43 %) 

0 ( 0 %) 

8 

21 ( 100 %) 

21 ( 100 %) 

14 ( 67 %) 

10 ( 48 %) 

1 ( 5 %) 

9 

49 ( 100 %) 

49 ( 100 %) 

30 ( 61 %) 

19 ( 39 %) 

1 ( 2 %) 

10 

165 ( 100 %) 

165 ( 100 %) 

110 ( 67 %) 

75 ( 45 %) 

1 ( 1 %) 

11 

552 ( 100 %) 

552 ( 100 %) 

355 ( 64 %) 

225 ( 41 %) 

10 ( 2 %) 

12 

2176 ( 100 %) 

2156 ( 99 %) 

1151 ( 53 %) 

727 ( 33 %) 

31 ( 1 %) 


Table 8. Distribution of I^{K) computations for HTW<i2. 
The entry in the cth row and reth column gives the number 
of knots K G HTWc for which it was necessary to compute 
the I^{K) in order to guarantee classification in HTW<i2. 


n\i 

2 

3 

4 

5 

6 

7 

3 

1 ( 100 %) 

1 ( 100 %) 

0 ( 0 %) 

0 ( 0 %) 

0 ( 0 %) 

0 ( 0 %) 

4 

1 ( 100 %) 

1 ( 100 %) 

1 ( 100 %) 

1 ( 100 %) 

0 ( 0 %) 

0 ( 0 %) 

5 

2 ( 100 %) 

2 ( 100 %) 

2 ( 100 %) 

2 ( 100 %) 

0 ( 0 %) 

0 ( 0 %) 

6 

3 ( 100 %) 

3 ( 100 %) 

3 ( 100 %) 

2 ( 67 %) 

0 ( 0 %) 

0 ( 0 %) 

7 

7 ( 100 %) 

7 ( 100 %) 

7 ( 100 %) 

3 ( 43 %) 

0 ( 0 %) 

0 ( 0 %) 

8 

21 ( 100 %) 

21 ( 100 %) 

16 ( 76 %) 

12 ( 57 %) 

1 ( 5 %) 

0 ( 0 %) 

9 

49 ( 100 %) 

49 ( 100 %) 

40 ( 82 %) 

24 ( 49 %) 

1 ( 2 %) 

0 ( 0 %) 

10 

165 ( 100 %) 

165 ( 100 %) 

142 ( 86 %) 

93 ( 56 %) 

6 ( 4 %) 

0 ( 0 %) 

11 

552 ( 100 %) 

552 ( 100 %) 

474 ( 86 %) 

301 ( 55 %) 

15 ( 3 %) 

0 ( 0 %) 

12 

2176 ( 100 %) 

2176 ( 100 %) 

1692 ( 78 %) 

1004 ( 46 %) 

48 ( 2 %) 

0 ( 0 %) 

13 

9988 ( 100 %) 

9972 ( 100 %) 

7198 ( 72 %) 

4426 ( 44 %) 

293 ( 3 %) 

2 ( 0 %) 


Table 9 . Distribution of I'^{K) computations for HTW<i3. 
The entry in the cth row and nth column gives the number 
of knots K G HTWc for which it was necessary to compute 
the I^{K) in order to guarantee classification in HTW<i3. 


algorithm | 5 . 3 | proposed in this paper computes a group presentation in a 
small fraction of time needed for the low-index subgroup computation, a 
reasonable strategy to minimize the total computation time is to test a few 
methods for each knot in the search of the possibly small number of critical 
cells. 

The problem is additionally complicated by the Tietze transformations 
used to simplify the group presentation. It turns out that the number of 
generators after transformations depends on some qualitative properties of 
the constructed discrete vector field. In Table [TT] we show the number of 
generators and relators before and after Tietze transformations in the case 
of four different strategies for the construction of the discrete vector field 
for knot K14,38437- Only in one case we get 3 generators and only in this 
case we are able to compute in several minutes. Our conclusion 
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n\i 

2 

3 

4 

5 

6 

7 

3 

1 (100%) 

1 (100%) 

1 (100%) 

1 (100%) 

1 (100%) 

0 (0%) 

4 

1 (100%) 

1 (100%) 

1 (100%) 

1 (100%) 

0 (0%) 

0 (0%) 

5 

2 (100%) 

2 (100%) 

2 (100%) 

2 (100%) 

0 (0%) 

0 (0%) 

6 

3 (100%) 

3 (100%) 

3 (100%) 

2 (67%) 

1 (33%) 

0 (0%) 

7 

7 (100%) 

7 (100%) 

7 (100%) 

5 (71%) 

2 (29%) 

0 (0%) 

8 

21 (100%) 

21 (100%) 

21 (100%) 

13 (62%) 

2 (10%) 

0 (0%) 

9 

49 (100%) 

49 (100%) 

46 (94%) 

26 (53%) 

3 (6%) 

1 (2%) 

10 

165 (100%) 

165 (100%) 

158 (96%) 

105 (64%) 

9 (5%) 

1 (1%) 

11 

552 (100%) 

552 (100%) 

523 (95%) 

329 (60%) 

26 (5%) 

0 (0%) 

12 

2176 (100%) 

2176 (100%) 

2001 (92%) 

1253 (58%) 

94 (4%) 

0 (0%) 

13 

9988 (100%) 

9988 (100%) 

8856 (89%) 

5494 (55%) 

451 (5%) 

3 (0%) 

14 

46972 (100%) 

46934 (100%) 

38092 (81%) 

23634 (50%) 

2231 (5%) 

21 (0%) 


Table 10 . Distribution of I'^{K) computations for 
HTW<i 4. The entry in the cth row and nth column gives 
the number of knots K G HTWc for which it was necessary 
to compute the I^{K) in order to guarantee classification in 
HTW<i4. 



Figure 2. Knot 1^14^38437- 

is that the number of generators does not determine the efficiency of the 
transformations. 

Nevertheless, for the calculations it is important to perform shaving and 
collapsibleSubset steps even for knots given as arc presentation (TH]. With¬ 
out the geometric simplifications we cannot get less than 4 generators for 
the knot in Figure The steps are also important for applications, where 
knots are placed in a big cubical grid, e.g. 3 D pictures of proteins. 
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Method 

Before Tietze 

After Tietze 

Reductions, cells order A 

g=4, r=4 

g=4, r=3 

Reductions, cells order B 

OQ 

II 

II 

g=3, r=2 

Coreductions, cells order A 

g= 6 , r =6 

g=4, r=3 

Coreductions, cells order B 

r—1 

II 

r—1 

II 

bjO 

CO 

II 

II 

b€ 
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